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ABSTRACT 

Using a complete, magnitude-limited sample of active galaxies from the Sloan Digital Sky 
Survey (SDSS) we show that the fraction of broad-line (Type 1) active galactic nuclei in¬ 
creases with luminosity of the isotropically-emitted [O III] narrow emission line. Our results 
are quantitatively in agreement with, and far less uncertain than, similar trends found from 
studies of X-ray and radio-selected active galaxies. While the correlation between broad-line 
fraction and luminosity is qualitatively consistent with the receding torus model, its slope is 
shallower and we therefore propose a modification to this model where the height of the torus 
increases slowly with AGN luminosity. We demonstrate that the faint-end slope of the AGN 
luminosity function steepens significantly when a correction for ‘missing’ Type 2 objects is 
made and that this can substantially affect the overall AGN luminosity density extrapolated 
from samples of more luminous objects. 
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1 INTRODUCTION 

Active galactic nuclei (AGNs) were once thought to be a rare phe¬ 
nomenon whose major cosmological impact was in their ability to 
be seen out to great distances. However, the discovery of super- 
massive black holes in the centres of nearby galaxies has led to the 
view that all galaxies went through an active phase. Furthermore, 
calculations have shown that the energy produced during this phase 
is sufficient to drive material out of the host galaxy (Silk & Rees 
1998; Fabian 1999). Feedback such as this can explain the tight¬ 
ness of the correlation between the masses of the black hole and 
host galaxy bulge (e.g., McLure & Dunlop 2002), and it has been 
suggested that such a process may perform a key role in shaping 
the galaxy luminosity function (Benson et al. 2003). 

Semi-analytic models of galaxy formation (e.g., GALFORM; 
Cole et al. 2000) are now being adapted to investigate the role of 
AGN-driven feedback, and determine whether a self-consistent pic¬ 
ture of both AGN and galaxy formation and evolution can be devel¬ 
oped. It is to be hoped that progress in this field may provide im¬ 
portant insights into some of the major outstanding issues in AGN 
research and, in particular, the question of AGN triggering. Early 
attempts at this (e.g., Kauffmann & Haehnelt 2000) have appeared 
encouraging in their ability to fit both the galaxy and QSO lumi¬ 
nosity functions but, by failing to make the distinction between the 
QSO and AGN luminosity functions, they fail to account for the 
possibly significant population of “Type 2 QSOs” which have been 
postulated to explain the shape of the Cosmic X-ray Background 
(CXB; e.g., Comastri et al. 1995). 

The cosmological importance of these objects is only ex¬ 
ceeded by their elusiveness, since their optical continua closely 
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resemble those of normal galaxies and therefore they cannot be 
identified by the colour-selection techniques used for QSOs (e.g., 
Schmidt & Green 1983). Radio selection has allowed the radio- 
loud population of Type 2 QSOs (i.e., radio galaxies) to be iden¬ 
tified, but radio sources comprise only a small fraction of the total 
AGN population. Hard X-ray selection has failed to locate the lumi¬ 
nous narrow-line AGN in the numbers predicted by model fits to the 
CXB but current X-ray telescopes are unable to find distant AGN 
with Compton-thick absorption, and such objects are known to ex¬ 
ist in abundance locally (e.g., Risaliti, Maiolino & Salvati 1999). 
Large spectroscopic surveys such as the Sloan Digital Sky Survey 
(SDSS; York et al. 2000) make it possible to identify AGN in an 
unbiased manner and therefore determine for the first time the frac¬ 
tion of Type 2 AGN. 

In the simple unified model for AGN, Type 1 and Type 2 ob¬ 
jects (i.e., those which display broad permitted lines, and those 
which show narrow lines only, respectively) differ only in terms 
of the angle which the observer’s line of sight makes with the axis 
of a dusty torus (see Antonucci 1993 for more details). If this angle 
is larger than some critical angle, 0c, the line of sight to the central 
regions (which includes the broad line region) is blocked by the 
torus and a Type 2 AGN is seen. QSOs are therefore only a subset 
of the AGN population and no cosmological distinction should be 
made. If 6c has no dependence on luminosity, then the QSO and 
AGN luminosity functions differ only in terms of their space den¬ 
sity normalization. 

However, the almost identical SEDs of QSOs just longward of 
1 pim (e.g., Kobayashi et al. 1993; Elvis et al. 1994) imply that the 
location of the inner wall of the torus is set by the distance from 
the nucleus at which dust reaches its sublimation temperature. This 
ides is supported by the observed time delays between variations 
in the optical and near-infrared continua (e.g., Eairall 9; Clavel, 
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Wamsteker & Glass 1989). In more luminous objects this distance 
is further away and, if the height of the torus remains constant, the 
critical angle 9o must increase, thereby leading to a luminosity de¬ 
pendence in the observed Type 1 AGN fraction. This is known as 
the ‘receding torus’ model and was first suggested by Lawrence 
(1991). 

Of course, the height of the torus need not be the same in all 
objects but, since the luminosity dependence of its inner radius is 
the result of basic physics, it would be even more surprising if its 
height varied in such a way as to produce a constant opening angle. 
A variation in the Type 1 AGN fraction with luminosity is there¬ 
fore to be expected, although the size of any such variation may be 
small. If it is large, however, then direct comparisons between the 
AGN luminosity function produced by models and the QSO lumi¬ 
nosity function derived from observations (such as that undertaken 
by Kauffmann & Haehnelt 2000) are seriously compromised. 

A number of authors have investigated possible evidence for 
the receding torus model, frequently using radio-selected AGN 
samples which are unbiased in terms of spectral type. These stud¬ 
ies (e.g.. Hill, Goodrich & DePoy 1996; Simpson, Rawlings & 
Lacy 1999; Willott et al. 2000; Simpson &. Rawlings 2000; Grimes, 
Rawlings & Willott 2004) have all supported the receding torus 
model, while Simpson (1998) showed how it could explain appar¬ 
ent differences in supposedly isotropic properties between radio 
galaxies and radio-loud quasars with the same radio luminosity. 
Recent deep hard X-ray surveys have enabled similar analyses to 
be performed on samples of radio-quiet AGN, and the same trend 
is seen (Ueda et al. 2003; Hasinger 2004). 

The case for the receding torus model appears strong, although 
most of the samples studies so far are relatively small. Alternative 
models for unification exist and include warped obscuring discs, 
broken shells of material, or an evolutionary scenario where the 
black hole is at first smothered but powers a wind which is even¬ 
tually able to blow away material and allow the QSO to be di¬ 
rectly observed. However, even if the receding torus model is in¬ 
correct, the number of Type 2 AGN need to be determined in or¬ 
der for AGN to fulfill their cosmological potential. In this paper, 
we construct separate [O III] emission line luminosity functions 
for broad and narrow-line AGN from the Second Data Release 
(DR2) of the SDSS which we use to investigate how the broad¬ 
line fraction varies with luminosity. The format of the paper is as 
follows. In Section 2 we describe how we have produced a sam¬ 
ple of emission-line galaxies from the SDSS DR2 catalogue, and in 
Section 3 we describe the steps taken to produce luminosity func¬ 
tions for both Type 1 and Type 2 AGN. In Section 4 we derive the 
Type I AGN fraction and discuss the implications of our findings. 
Finally in Section 5 we present a summary. Throughout this pa¬ 
per we adopt a flat cosmology with flm = 1 — Da = 0.3 and 
Ho — 70kms“^ Mpc“^. 


2 SAMPLE SELECTION 

The SDSS DR2 spectroscopic catalogue covers 2627 deg^ down 
to a Petrosian magnitude of r' < 17.77 (in addition to colour- 
selected QSO candidates which we exclude by means of this mag¬ 
nitude cut). Spectra are taken through 3-arcsecond diameter fibres, 
and so we must impose distance (i.e., redshift) limits on our sam¬ 
ple. The fibres will not encompass the entire emission line region 
for nearby AGN, since this can be ~ 1 kpc in size (e.g., Wilson & 
Tsvetanov 1994), while the dilution by starlight in distant galaxies 
will make the detection of emission lines more difficult. Kauffmann 


et al. (2004) show that a Seyfert 2 sample with a luminosity limit 
^[oiii] > 10^ 1/0 and redshift limits 0.02 < a < 0.30 does not 
suffer from aperture biases, and we therefore adopt the same con¬ 
straints for our sample. This produces a catalogue of 4304 objects. 

We do not apply any reddening corrections to the line lumi¬ 
nosities. This is partly for practical reasons, since there is no reli¬ 
able way to determine the reddening to the narrow-line region in 
broad-line objects. This is always a problem for studies which aim 
to compare the emission line properties of Type 1 and Type 2 AGN, 
although the fact that [O III] emission is seen to be isotropic when 
no reddening correction is applied (Mulchaey et al. 1994) suggests 
that differences between the mean reddening to the two types of 
galaxy must be small. Dahari & De Robertis (1988) estimate that 
the median difference is E{B — H) « 0.2, corresponding to a fac¬ 
tor of ~ 1.7 at 5007 A. We will address whether this could have an 
impact on our results later. 

We also do not correct our emission line fluxes for underlying 
stellar absorption. Since all our galaxies have very luminous [O III] 
emission lines, they will also have luminous H/7 emission which 
is expected to dominate over any absorption. Wake et al. (2004) 
find that the observed AGN fraction is not sensitive to whether this 
correction is made. 

3 ANALYSIS 

3.1 Spectral classification 

The SDSS spectroscopic pipeline produces a classification for all 
objects, based on principal component analysis (Schlegel et al., in 
preparation). For our sample of objects, this classification is either 
“QSO” (890 objects) or “galaxy” (3414 objects), which might be 
expected to correspond to AGN Types 1 and 2, respectively. How¬ 
ever, while all objects classified as QSOs are expected to display 
broad emission lines, the galaxy classification will include star- 
burst galaxies as well as objects classified as Type 1.8 or 1.9 (where 
broad Ha emission is seen, but not broad H/3) and possibly some 
low-luminosity Seyfert Is. Since broad emission components will 
preclude the use of a line-ratio diagram to discriminate between 
starburst and AGN galaxies (Baldwin, Phillips & Terlevich 1981), 
we elect to classify all 4304 galaxies by eye. 

It is essential that the classification we choose is not depen¬ 
dent on the emission line luminosity. Since the broad and narrow 
emission line luminosities are strongly correlated, this requires that 
reddened objects where the broad emission component is near the 
limit of detection be included with the narrow-line objects. Objects 
which display broad wings on the Ha line (‘Type Lx’ galaxies) are 
therefore grouped with the Type 2 galaxies, while those which have 
clear broad Ha and H/3 emission lines are classified as Type Is. 

One object was discovered to have a large noise spike at the 
expected location of [O III] A5007 and was removed from the sam¬ 
ple, since the low luminosities of the visible emission lines and 
the absence of the A4959 line indicate that its true line luminosity 
would be far below our lower limit. 

We find that 63 objects classified as QSOs (7.1 per cent) are 
not QSOs at all, but display only narrow emission lines and ap¬ 
pear to be Type 2 galaxies. We therefore reclassify them as such. 
We also find that 63 objects classified as galaxies (1.8 per cent) are 
genuinely QSOs, with clear broad line emission from both Ha and 
H/3, and we reclassify these objects. A further 102 (3.0per cent) 
show quite pronounced broad wings on the Ha line only, with a 
similar number showing weaker broad wings whose detectability 
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Figure 1. Emission-line ratio diagram for narrow-line objects (dots) and 
Type 1.x galaxies (filled circles). The a and /3 axes represent the ridge lines 
from the intersection of the AGN and starburst loci. The dotted and dashed 
lines are the classification criteria suggested by Kewley et al. (2001) and 
Kauffmann et al. (2004), respectively. 


would be a function of the signal-to-noise ratio in the spectrum. 
A subsample of 100 objects was classified by Dr C. Done, and it 
was only for objects in this last category where disagreement was 
found. The total fraction of SDSS-classified ‘galaxies’ which we 
believe display broad line emission is therefore consistent with the 
figure of 8 per cent quoted by Kauffmann et al. (2004). 

While the reason for the misclassification of Type 2 galaxies 
as QSOs is unclear, the Type Is which are classified as galaxies 
typically have either weak broad line emission, or have very nar¬ 
row permitted lines (~ 500kms“^). Our interpretation of these 
as bona fide broad-line objects is based on the presence of Fe II 
pseudo-continuum emission between Ft 7 and H/3, the lack of a 
strong 4000-A break or Balmer jump, and a continually rising spec¬ 
trum at short wavelengths. 


3.2 Separation of Seyfert and starburst galaxies 

We wish to construct comparable Type 1 and Type 2 samples, so 
that their luminosity functions (LFs) are properly representative of 
the population as a whole. Unfortunately, there is no straightfor¬ 
ward way to do this. The objects we have classified as Type Is all 
contain an AGN component, since the broad emission lines we ob¬ 
serve cannot be produced by stellar processes, but they may also 
possess a significant contribution to their emission line luminosi¬ 
ties from star formation activity. Objects where the line emission is 
almost entirely powered by star formation may not be included in 
this sample, however, even if we have a direct view of the broad line 
region, since the broad line emission will be weak. We therefore 
need to remove pure starburst galaxies from the sample of narrow¬ 
line objects, but also objects where the AGN contribution is so low 
that the objects would not have been classified as Type Is if we had 
been able to see the BLR. 

Kauffmann et al. (2004) use an emission-line ratio diagram 



a 

Figure 2. FWHM of the [O III] emission lines plotted against distance along 
the a axis of Fig. Q for narrow-line objects (dots) and Type 1.x galaxies 
(filled circles). 


of the type pioneered by Baldwin et al. (1981) to separate Seyferts 
from starburst galaxies, and we propose a similar method. Since the 
appropriate diagram makes use of the narrow Balmer line fluxes, it 
cannot be used for Type 1 AGN, where these data are unavailable. 
By plotting the Type 1.x galaxies on the same diagram, however, 
we can estimate what regions the Type 1 population would occupy. 
Fig.Qdisplays these data. 

Since our ‘Type 1.x’ classification requires that no broad H/3 
emission is seen, the location of these points along the ordinate 
should be representative of the Type 1 population. Many of the 
Type 1.x galaxies lie in the same region as the Type 2s, because 
the broad Fla emission is weak and the SDSS measurement there¬ 
fore only includes the narrow-line flux. However, the tendency for 
these some of objects to have lower [N Il]/Ha ratios than the Type 2 
population is clear, as the Ha flux measurement for these includes 
the broad line emission. What is also clear is that some galaxies 
containing an AGN have low [O IIl]/H/3 ratios, and therefore any 
cut which excludes starburst galaxies will also exclude some bona 
fide AGN. 

We attempt to find the appropriate criteria for separating 
Type 2 galaxies from pure starbursts by defining two (non- 
orthogonal) axes as shown in Fig.Q The axes meet where the loci 
of the starburst and AGN sequences meet. The ‘a’ axis follows the 
ridge line of the AGN sequence, while the ‘/3’ axis follows the ridge 
line of the starburst sequence. In this system, the demarcation cri¬ 
teria proposed by Kewley et al. (2001) and Kauffmann et al. (2004) 
correspond to a « 0.5 and a « 0, respectively. 

We also consider the widths of the narrow emission lines 
(specifically, [O III] A5007) to help the differentiation between 
Seyfert and starburst galaxies, since it is known that the narrow 
lines are broader in Seyfert galaxies than in starbursts (Shuder & 
Osterbrock 1981). This is demonstrated to be the case in Fig. 
where we have plotted the FWHM of the [O III] line against a as 
determined from Fig.Q For the narrow-line galaxies, we have taken 
the value of a from their location in the (a,/3) coordinate system. 
Since the [N Il]/Ha ratio for the Type 1.x galaxies is affected by 
broad line emission, we compute values of a for these galaxies by 
determining where a galaxy with the same [O IIl]/H/3 ratio would 
lie on the a-axis. 

The overall distributions of the Type 2 and Type 1.x galaxies 
are inconsistent at the ^ 99.99% level, using the two-dimensional 
Kolmogorov-Smirnov test of Peacock (1983). However, Fig. 
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clearly shows two populations of narrow-line galaxies - a popula¬ 
tion with broad (FWHM > 250 kms“^) emission lines extending 
to large values of a which appears to be well traced by the Type 1 .x 
galaxies; and a population of galaxies with narrower emission lines 
and a « 0. These two populations can be associated with Seyfert 
and starburst galaxies, respectively, and it is necessary to find some 
way to discriminate between the two classes. 

Since the two distributions obviously differ greatly in the num¬ 
ber of galaxies in the lower left comer of Fig.|^ we investigate dis¬ 
crimination criteria which remove these objects from the narrow¬ 
line sample. We construct a family of four-parameter selection cri¬ 
teria which require that a and the emission line FWHM are both 
greater than certain values, and also that the points satisfy the rela¬ 
tion log FWHM > m{ao — a) where m > 0. We then compare 
the samples of narrow-line galaxies selected by each of these cri¬ 
teria with the complete sample of Type 1 .x galaxies in the param¬ 
eter space of Fig. again using a two-dimensional K-S test. We 
find that the two distributions can be made consistent at better than 
90per cent confidence by making a simple cut with a > 0.20. No 
improvement is made by imposing an additional cut based on line 
width, and only a very marginal gain results from adding a further 
selection criterion which is a function of both parameters. This se¬ 
lection lies approximately halfway between the criteria suggested 
by Kewley et al. (2001) and Kauffmann et al. (2004), as expected. 
The implication is that the region between their two sets of criteria 
represents a mixing line, with our limit being the point at which 
the AGN contribution becomes significant (i.e., when it is visible 
in the spectrum even if lightly-reddened). In the sample of 4303 
emission-line galaxies there are 851 Type 1 AGN, 141 Type 1.x 
AGN, 2184 Type 2 AGN, and 1127 starburst galaxies (using the a 
criterion to classify the latter two classes). 

3.3 Comparison of luminosity functions 

Since the SDSS DR2 spectroscopic survey is a complete 
magnitude-limited survey. Type 1 galaxies will appear in relatively 
greater numbers because of the additional flux from the nuclear 
non-stellar continuum. We can however account for this by con¬ 
structing LFs using the 1/14 estimator (Schmidt 1968; Rowan- 
Robinson 1968). We then compare the Type 1 LF with the LFs of 
three subsamples of Type 2 galaxies, constructed by applying dif¬ 
ferent selection criteria to the entire sample in an attempt to exclude 
the starburst galaxies. Two of these subsamples, named Type 2-Ka 
and Type 2-Ke, are produced by applying the line ratio selection 
criteria of Kauffmann et al. (2004) and Kewley et al. (2001), re¬ 
spectively, while the third, named Type 2-a, is constructed using 
the a > 0.2 criterion derived above. 

The three Type 2 LFs are compared with the Type 1 LF in 
Fig. 13 It is clear that, although the numbers of each type of Seyfert 
are similar at L[oiii] 10^^'® W, Type 2s outnumber the Type Is 
at lower luminosities by a factor of ~ 5. Although there are dif¬ 
ferences between the numbers of lower luminosity Type 2s in the 
three samples, these differences are small compared to the differ¬ 
ence between these samples and the number of Type Is. 

The different shapes of the LFs cannot be explained by a dif¬ 
ference in the mean reddening of the Type 1 and Type 2 samples, 
since a reddening correction is simply a horizontal shift. The differ¬ 
ent shapes also cannot be due to problems with the Seyfert/starburst 
separation. The criterion defined by Kewley et al. (2001) is con¬ 
servative in the sense that it should exclude all starburst galaxies, 
so the excess of Type 2 galaxies over Type 1 galaxies at lower 
[O III] luminosities must be genuine, since even the Type 2-Ke 



Figure 3. Luminosity functions for Type 1 galaxies (filled squares) and 
samples of Type 2 galaxies selected from the nan'ow-line objects using the 
criteria of Section 3.2 (open squares), Kewley et al. (2001; open circles), 
and Kauffmann et al. (2004; open triangles). Also shown (filled triangles) 
is the luminosity function of the Type 1.x galaxies, although these are also 
included in all the Type 2 samples. 

sample exceeds the Type 1 sample. Conversely, the difference also 
cannot be due to an artificial deficit of high luminosity Type 2 
galaxies, since all three samples have effectively the same LFs at 

i[oiiii>10®"®W. 

A final possibility is perhaps that a significant fraction of the 
lower luminosity Type 1 galaxies have been misclassified, pre¬ 
sumably because the broad lines are weak and below the detec¬ 
tion threshold. This seems unlikely since the broad lines are much 
stronger than the narrow lines. In our Type 1 sample, the median 
[O IIl]/Ha ratio is 0.15, while 90 per cent of objects have a ratio 
less than 0.5. We therefore expect all but a small fraction of Type 1 
galaxies to have readily identifiable broad Ha emission (although 
the [O III] lines in our sample only need to be detected at S/N > 5, 
most are detected at much higher significance - all but 20 are de¬ 
tected at S/N > 10, and over half have S/N > 20). Furthermore, 
we also display the LF of the Type Lx galaxies in Fig.|3 and there 
is clearly no bias against lower luminosity objects. We therefore 
conclude that the luminosity functions of Type 1 and Type 2 galax¬ 
ies are genuinely different. 

For the remainder of this paper, we use only one Type 2 lumi¬ 
nosity function. We adopt the Type 2-a sample, and add in quadra¬ 
ture to the Poisson uncertainty the standard deviation of the three 
samples. 

4 DISCUSSION 
4.1 The Type I fraction 

The different luminosity functions for our samples of Type 1 and 
Type 2 galaxies can also be described as a luminosity dependence 
of the Type 1 AGN fraction. In the context of the receding torus 
model, such a change is described by 

/i =l-[l + 3L/Lo]-°® (1) 

(e.g., Simpson 1998), where Lq is the luminosity at which the num¬ 
bers of Type 1 and Type 2 AGN are equal (i.e., an opening angle 
of = 60°). The ‘luminosity’ of concern here is the ultraviolet- 
optical radiation which heats the dust, which is difficult to measure 
directly, but the constant equivalent width of the [O III] emission 
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Figure 4. Type 1 AGN fraction as a function of [O III] luminosity, from this 
paper (solid squares), Ueda et al. (2003; open squares), Hasinger (2004; 
open circles), and Grimes et al. (2004; open triangles). The data from Ueda 
et al. and Hasinger have been converted to L[o III] described in the text. 
The solid line shows the best-fitting ‘standard’ receding torus model (Equa- 
tionQ while the dashed line shows the best-fitting model where the torus 
height is allowed to vary (Equation|^, and the dotted line shows the best¬ 
fitting model where the proportionality between L^q jjjj and Lpad is bro¬ 
ken (Equation^. 

line in QSOs (e.g. Miller et al. 1992; see also Simpson 1998) sug¬ 
gests that it should be a good tracer. Since the receding torus model 
is able to explain a wide range of observed AGN properties (Simp¬ 
son 2003 and references therein; see also Grimes et al. 2004), we 
investigate whether it is also able to explain the results of the pre¬ 
vious section. 

We complement our data from similar analyses in the litera¬ 
ture, which typically extend to higher luminosities. Both Ueda et al. 
(2003) and Hasinger (2004) measure the Type 1 AGN fraction in 
hard X-ray selected samples, which should be relatively free from 
orientation biases except for Compton-thick AGN, whose existence 
will case the Type 1 fraction to be systematically overestimated. 
About half of nearby low luminosity Type 2 AGN are Compton- 
thick (Risaliti et al. 1999) and, if this fraction is true at all lumi¬ 
nosities then the Type 1 fraction will be overestimated by about 
50 per cent. We do not make any correction for Compton-thick ob¬ 
jects but note that the Type 1 fractions derived from X-ray-selected 
samples should be considered as upper limits to the true fractions. 
We convert the 2-10 keV luminosity (corrected for intrinsic absorp¬ 
tion) used by these authors into [O III] luminosity using the mean 
ratio for Seyfert galaxies from Mulchaey et al. (1994), i.e., 

i[0 III] = 0.015 X I/2-lOkeV • (2) 

We also add the data of Grimes et al. (2004) from the 3CRR, 6CE, 
and 7C flux-limited samples of radio sources. The results are shown 
in Fig. |4| where a clear and highly significant (> 99.99 per cent 
confidence) correlation between the Type 1 fraction and luminos¬ 
ity is seen for which the receding torus model provides an attrac¬ 
tive, but not unique, explanation for this correlation. We note that 
the points from Grimes et al. (2004) tend to be slightly lower than 
those from the other works and this might reflect a difference in the 
torus properties of radio-loud and radio-quiet AGN, although the 
discrepancy is not highly significant. 

Although the four samples plotted in Fig.^all show the same 
trend, we concentrate for now on the data from our analysis alone 
since there are less well understood selection effects in the X-ray 


samples (see Treister et al. 2004). These cover more than two orders 
of magnitude in L[o m] but are not well-fit by the standard receding 
torus model, with the best fit having = 65.8/9 for Lq — 

1034.90 

We therefore suggest a modification to the standard receding 
torus model, where the height of the torus depends on the luminos¬ 
ity of the AGN. Such a dependence could plausibly be modulated 
via the mass of the black hole and the strength of its associated 
gravitational potential. We parametrize this dependence as oc 
and EquationQtherefore becomes 

/i = 1-[l + 3(L/Lo)i-"«]-° ". (3) 

This produces a more acceptable fit (of jv = 12.1/8) for 
I/O = 10^®'®^ W and ^ = 0.23, which we also show in Fig.|3] 
Arshakian (2004) made a similar parametrization when studying 
the luminosity-dependent variations of QSO fraction and projected 
linear size in a radio source sample. He defined a parameter A 
such that tan 0c oc (therefore A = 0.5 — 0 and found 
As = 0.35 ± 0.09 and A{ = 0.26 ± 0.02 when considering the 
linear sizes and QSO fractions, respectively. The latter value is in 
excellent agreement with the best-fitting value of / we determine 
here. 

An alternative explanation for this modification is that I[o m] 
is not a good substitute for Irad. the radiative luminosity which 
heats the dust in the torus. If I[o m] oc then the stan¬ 

dard receding torus model would fit the data. However, such an 
extreme disproportionality can be excluded based on a number of 
arguments such as the near-constancy of [O III] equivalent width 
in QSOs, (Miller et al. 1992), the proportionalities between [O III] 
luminosity and other quantities (e.g., Rawlings & Saunders 1991; 
Mulchaey et al. 1994; Grimes et al. 2004), and the dependence of 
[O III] emission on ionization parameter (Simpson 1998). 

A more subtle variation in the relationship between lead and 
I[o III] may provide an explanation. In the standard unified model 
for AGN, only a fraction ^ fi of the ionizing radiation can escape 
to illuminate the narrow-line region on scales > 1 pc. This explains 
the (bi)conical emission-line structures seen in some Seyfert galax¬ 
ies (e.g. Pogge 1989; Haniff, Wilson & Ward 1991; Simpson et al. 
1997), and suggests that I[o m] Irad/i- In such a scenario, 

/i = 1 — [1 + l-5I[o iii]/(Io/i)l , (4) 

which can be rearranged to give a cubic equation in /i whose roots 
are functions of I[o m]- This model is also shown in Fig. 0 and 
provides an acceptable fit (i.e. cannot be rejected at 95 per cent con¬ 
fidence) to our data {yf /v = 16.0/9 for Lq = W). 

Fig. 131 suggests that we should favour the torus described by 
Eqn|3|over that from Eqn|4|by virtue of the quality of fit to the 
X-ray datapoints. However, we caution against leaping to such a 
conclusion and suggest that some of the most luminous QSOs may 
have been misclassified as Type 2 objects. While this seems coun¬ 
terintuitive, such objects will inevitably be at the highest redshifts 
and therefore the optical spectroscopy performed on them will sam¬ 
ple the rest-frame ultraviolet. Even relatively small amounts of dust 
could be sufficient to extinguish the broad emission lines and give 
such objects the appearance of a Type 2 AGN in the UV. While 
there should be very few luminous AGN which are lightly reddened 
because the observer’s line of sight just grazes the torus (see Hill 
et al. 1996), the required extinctions may arise from dust in the 
host galaxy which is not associated with the active nucleus. The re¬ 
peated reclassifications of ‘Type 2 QSOs’ after near-infrared spec¬ 
troscopy reveals broad Ha emission (e.g., Halpem, Eracleous & 
Forster 1998; Halpem, Turner & George 1999; Akiyama, Ueda & 
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Ohta 2002) suggests that there may be some truth in this. On the 
other hand, the QSO fraction may be overestimated if there exists 
a significant population of luminous AGN with Compton-thick tori 
which are not detected in hard X-ray surveys. What is required is 
a systematic near-infrared study of either spectroscopy of imaging, 
similar to that which Simpson et al. (1999) performed for a ~ 1 
powerful radio sources. This should be possible in the very near fu¬ 
ture with the combination of Chandra or XMM-Newton and Spitzer 
imaging in a number of fields. 

Having suggested that both modifications to the receding torus 
model fit the data equally well, we return to the issue of whether 
breaking the proportionality between the ionizing and [O III] 
emission-line luminosities is in conflict with observations. First, 
this does not affect the isotropy of the [O III] emission, since the 
emission-line luminosities of Type 1 and Type 2 AGN with similar 
ionizing luminosities should still be the same. However, it predicts 
that the observed equivalent width of [O III] should be proportional 
to fi (i.e., increase with increasing AGN luminosity), whereas ob¬ 
servations show it to be constant with luminosity (Miller et al. 
1992). However, over the range of luminosities studied by Miller et 
al. (approximately two orders of magnitude), /i is predicted to vary 
by a factor of about 2, which is equal to the la scatter in the rest- 
frame equivalent width. Furthermore, this effect will be lessened 
by those objects seen at small polar angles where the continuum, 
but not the emission-line, luminosity is enhanced. A quantitative 
analysis of this effect depends on the assumed angular and wave¬ 
length dependence of the ionizing continuum and is therefore rather 
model-dependent and beyond the scope of this paper. However, we 
do predict that the scatter in the equivalent width of the [O III] line 
should increase with luminosity. 


4.2 Implications for AGN evolution 


Studies of AGN evolution have so far primarily relied on samples of 
QSOs since these objects are readily visible out to large distances. 
Obviously such samples contain only a fraction of the AGN at any 
redshift or luminosity since Type 2 objects do not meet the colour 
and/or morphological selection criteria (while hard X-ray surveys 
are less biased against Type 2 AGN, they presently lack the com¬ 
bination of depth and area coverage to properly sample the AGN 
luminosity function in an unbiased manner). We therefore wish to 
determine how the luminosity-dependent QSO fraction might af¬ 
fect the conclusions drawn from studies of QSO evolution. 

Since there is a significant stellar component to the broad-band 
magnitudes of our Type 1 galaxies which is not the case for ob¬ 
served samples of QSOs, we estimate the non-stellar continuum 
luminosities from our [O III] luminosities.We assume a mean rest- 
frame equivalent width (relative to the non-stellar continuum) for 
the A5007 line of e = 24 A (Miller et al. 1992), and a spectral shape 
Si, QC (Vanden Berk et al. 2001; the result is extremely in¬ 

sensitive to the shape of the spectrum due to the small wavelength 
difference between the [O III] line and the B band) in the rest-frame 
optical. It then follows that 

Mb = -22.0 - 2.51og(L[oiii]/10®® W) + 2.51og(£/24A) (5) 


and Mhj = Mb — 0.09 for the canonical QSO spectrum (Blair & 
Gilmore 1982). 

Our Type 1 LF can be well-fit with the two power-law 
parametrization of the 2QZ team. 




2<f>(L*) 

(L/L*)-“ -F (L/L*)-/3 


( 6 ) 



Figure 5. Our Seyfert (Types 1 and 2) luminosity function converted to 
Mhj (crosses) and plotted with the modified QSO optical luminosity func¬ 
tions of Groom et al. (2004; all other symbols). A two power-law model 
with pure luminosity evolution has been fitted to the data, as described in 
the text, and the results are shown by the solid curves. 


(e.g.. Groom et al. 2004; note that the normalization factor of 2 is 
missing from papers by this team). A good fit lod jv — 4.8/6) 
is found with the parameters a = —2.86, /3 = —1.06, L* = 
3.5 X 10 ^* W, and = 2.24 x 10“’’ Mpc“® dex“’. The power 
law indices are similar to those derived by Groom et al. (2004), who 
found a = —3.25 and j3 = —1.01. However, a better comparison 
would be to the combined (i.e.. Types 1 and 2) Seyfert LF since, 
while the Groom et al. (2004) LFs are relatively insensitive to the 
lack of Type 2 objects, our Type 1 LF is a poor representation of 
the overall AGN LF. The fit to the combined Seyfert luminosity 
function is improved jv = 3.3/6), primarily due to the larger 
error bars resulting from the uncertain classification of narrow-line 
objects, with a = —3.11, (5 = —1.59, L* = 4.4 x lO®"’W, and 
= 5.25 X 10“'^ Mpc“^ dex“A 

There is moderately good agreement between the values of 
M* and from our combined Seyfert sample and the extrap¬ 
olation of the exponential pure luminosity evolution of Groom 
et al. (2004). We obtain = -21.20 and = 1.31 x 

10“® Mpc“® mag“’ compared to = —21.46 and d?* = 
1.84 X 10“® Mpc“® mag“’ from their work (the agreement is 
less good for the polynomial evolution, where a value of = 
—22.15 is predicted for our sample). The lower space density we 
derive cannot be due to incompleteness in our sample of Seyfert 
galaxies since our spectroscopic selection will include all such ob¬ 
jects (accounting for the uncertainty in the starburst-Seyfert classi¬ 
fication), whereas the optical colour selection used to select QSOs 
will exclude a population of red QSOs which may exist in signif¬ 
icant numbers (Webster et al. 1995; Benn et al. 1998; Richards et 
al. 2003). 

We also try to simultaneously fit our Seyfert LF and the six 
QSO LFs of Groom et al. (2004), after using the modified reced¬ 
ing torus model to account for the missing Type 2 objects. We use 
their two power-law model with an exponential form for the lumi¬ 
nosity evolution, although we fit the value of M* for our Seyfert 
LF as a separate parameter since an extremely poor fit is obtained 
otherwise. Following Groom et al., we use only the Mbj < —22.5 
bins, and we show the best fit in Fig.|3 While this fit is unaccept¬ 
ably poor in a formal sense jv — 175/58), it does reproduce 
the gross characteristics of the evolving LF and, in particular, the 









The luminosity dependence of the Type 1 AGN fraction 1 


power law indices differ from those derived from the Type 1 (QSO) 
LF alone. 

Unsurprisingly, the value of a. is least affected. From our ear¬ 
lier calculation that Lq = W, we conclude that QSOs out¬ 

number Type 2 objects where Ms < —23.3 and therefore the bright 
end of the luminosity function is not strongly affected by the ab¬ 
sence of Type 2 objects and the power law slope is only moderately 
steepened by their inclusion - we find a. = —3.42 as opposed to 
a = —3.25. However, the faint end slope is made much steeper 
(j3 = —1.41 compared to /3 = —1.01), and is closer to the faint- 
end slope of the hard X-ray luminosity function {j3 = —1.86; Ueda 
et al. 2003). This steepening has a large effect on the contribution 
from faint AGN to the overall AGN luminosity. Integrating the two 
power-law model, the total luminosity from AGN with L > L* 
is ~ QA^*L*. The total luminosity from AGN fainter than L* 
does not converge for /3 < —1 and so is dependent on the low- 
luminosity cutoff, but in the decade below L*, it is ~ 2.34>*L* 
for j3 ~ —1.0 and nearly twice as much for /3 « —1.5. Wolf et 
al. (2003) find a downturn in the shape of the LF at faint abso¬ 
lute magnitudes (Mb > —22, not probed by 2QZ) which allows 
the luminosity density integral to converge. However, the increased 
number of Type 2 AGN predicted at these luminosities significantly 
affects the strength of this downturn and therefore the contribution 
by faint AGN to the overall luminosity density. Such objects would 
be fainter than the low-luminosity Type 1 AGN studied by Wolf et 
al. (2003) and therefore could not have yet been spectroscopically 
identified in deep X-ray surveys such as Barger et al. (2003). Given 
that a putative further population of Compton-thick Type 2 AGN 
would not even appear in such X-ray surveys, there is a distinct 
likelihood that the importance of low luminosity AGN has been 
significantly underestimated. 

Finally, while attempting to tie the present-day black hole 
mass density to the evolution of QSOs, Yu & Tremaine (2002) note 
that their results are dependent on the prevalence of obscured low- 
luminosity AGN. In particular, the accretion efficiency can be made 
roughly independent of AGN luminosity if the obscured AGN frac¬ 
tion is larger for low mass black holes, and hence less luminous 
AGN, as predicted by the receding torus model. 

4.3 The Cosmic X-ray Background 

While the predictions of the receding torus model agree well with 
the observed number of X-ray-selected Type 2 AGN, models for 
the CXB (e.g., Comastri et al. 1995; Wilman, Fabian & Nulsen 
2000) require the presence of an additional Compton-thick popu¬ 
lation which is not detectable in such surveys. As discussed previ¬ 
ously, this means the true Type 1 fraction is lower than the observa¬ 
tions suggests, and hence there is an apparent discrepancy between 
the receding torus model and models for the CXB. 

However, the number of Compton-thick AGN proposed is a 
relatively small fraction of the total AGN population. Ueda et al. 
(2003) fit the CXB by proposing that there are as many AGN 
with lO^'^cm”® < A^h < 10^® cm“® as with 10^® cm“® < 
Nh < 10^'*cm“®, but there is a strong luminosity dependence 
in the shape of the A^h distribution. Therefore, while there are 
many Compton-thick AGN at low redshift/luminosity (Risaliti et 
al. 1999), there are many fewer at the redshift/luminosity which 
produces the bulk of the CXB. 

In addition. Type I AGN with significant X-ray obscura¬ 
tion are known to exist, the original member of this class being 
MR 2251 — 178 (Halpern 1984; see also Page et al. 2004), so the 
Type 1 fraction need not be the same as the fraction of X-ray unab¬ 


sorbed AGN. Although the column densities of these objects tend 
not to exceed ~ 10^^ ® cm“^, Willott et al. (2003, 2004) have 
found objects with greater absorption (Wh > 10^® cm^) whose 
optical/infrared spectra are consistent with lightly-reddened QSOs. 
It is therefore possible that the most heavily absorbed AGN also 
suffer significant non-torus reddening and may have escaped de¬ 
tection in both X-ray and optical QSO surveys. 

Finally, a large discrepancy between the receding torus model 
and the observations of hard X-ray-selected AGN occurs only at 
very high luminosities (I/[oiii] > 10®® W, or Lx > 10®® W) and 
this population contributes almost negligibly to the overall CXB 
(e.g., figs 16 and 17 of Ueda et al. 2003). We conclude therefore 
that the receding torus model does not conflict with our current 
understanding of the Cosmic X-ray Background. 


5 SUMMARY 

We have constructed separate luminosity functions for Seyfert 1 
and Seyfert 2 galaxies, using the Second Data Release of the Sloan 
Digital Sky Survey. We have shown that these luminosity functions 
have different shapes, resulting in an increase in the Type 1 AGN 
fraction with luminosity which is in quantitative agreement with 
that determined from hard X-ray-selected samples. We show that 
the form of this luminosity dependence is inconsistent with the 
standard receding torus model, and propose a modification where 
the height of the torus varies with luminosity as /i oc L®'^®, al¬ 
though we note that a more robust spectral classification of dis¬ 
tant luminous AGN is desired. Assuming that this modification to 
the receding torus model provides a reasonable correction factor 
to account for the ‘missing’ Type 2 AGN in QSO surveys, we have 
shown that the contribution to the cosmic accretion luminosity from 
faint AGN may have been significantly underestimated. 
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